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基于 粒子 群 算法 的 三 维 CAD 模型 相似 性 计算 ， 
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摘 要 : 为 了 更 加 准确 地 度量 两 个 模型 之 间 的 形状 差异 ， 提 出 了 一 种 基于 粒子 群 的 模型 相似 性 计算 方法 。 利 用 面 的 组 
成 边 数 来 构造 面相 似 度 和 矩阵， 通过 粒子 群 算法 对 该 矩阵 进行 搜索 ， 得 到 了 两 个 模型 之 间 的 最 优 面 匹 配 序列 。 根 据 这 个 
最 优 面 匹配 序列 ,从 面相 似 度 矩阵 中 提取 对 应 的 面相 似 性 值 。 通 过 累积 面 之 间 的 相似 性 来 计算 模型 之 间 的 整体 相似 性 

以 此 为 基础 来 度量 模型 之 间 的 差异 。 实 验 结果 表明 : 该 方法 能 够 准确 地 度量 两 个 模型 之 间 的 相似 程度 
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Similarity calculation of 3D CAD model based on particle swarm algorithm 


Gao Xueyao’, Chen Yunan’, Zhang Chunxiangt 
(a. School of Computer Science & Technology, b. School of Software, Harbin University of Science & Technology, Harbin 
150080, China) 


Abstract: In order to measure shape difference of two models accurately, this paper proposed a method for calculating model 
similarity based on particle swarm. It used edge number of face to construct face similarity matrix, and used particle swarm 
algorithm to search this matrix to get an optimal face matching sequence between two models. It extracted correspondent face 
similarity values according to this face matching sequence from face similarity matrix. This paper accumulated face similarities 
to calculate two models' similarity and measured shape difference of two models. Experimental results show that this method 
can evaluate similarity between two models accurately. 


Key Words: shape difference; model similarity; particle swarm algorithm; matching sequence 


| 


0 引言 


非 线性 凝聚 层次 聚 类 算法 来 实现 属性 图 聚 类 。 同 时 ， 引 入 增 量 
模型 的 动态 归 类 方法 ， 来 实现 三 维 CAD 模型 可 重用 区 域 的 

在 三 维 CAD 模型 检索 中 ， 模 型 相似 性 计算 起 着 重要 的 作 ”，” 效 聚 类 辐 。Qin 利用 分 层 特 征 本 体 和 本 体 映 射 , 来 生成 语义 描述 
用 。 它 关系 到 模型 检索 的 效率 和 可 靠 性 。 刘 志 通 过 获取 三 维 模 。 符 。 以 语义 描述 符 为 基础 进行 本 体 推理 ， 使 检索 系统 获得 更 好 
型 较 优 视点 集 ， 来 这 染 三 维 模型 的 混合 轮廓 线 视 图 。 利 用 提取 ”的 性 能 外 。Kim 给 出 了 一 种 基于 自由 草图 的 建 模 方法 。 利 用 形 
的 Gabor 边缘 响应 特征 ， 去 构建 特征 库 。 从 特征 库 中 ， 使 用 视 。 ” 状 分 布 来 比较 多 分 辨 率 模型 ， 使 用 简单 的 形状 查询 模型 来 提高 
觉 词 袋 方法 来 检索 相似 的 模型 趾 。 范 壮 采 用 多 层次 局 部 视图 对 检索 精度 由。Chen 给 出 了 三 个 组 合 特 征 描述 符 和 一 个 基于 类 的 
树木 模型 进行 组 合 描述 ， 根 据 视图 所 对 应 的 树木 语义 信息 ,来 ” ”特征 描述 符 , 并 将 其 应 用 于 模型 检索 过 程 中 四 。Tao 把 面 邻接 图 
实现 视图 特征 的 转换 , 从 而 完成 由 粗略 到 精细 的 多 阶段 检索 中。 分割 为 目 、 止 和 平面 区 域 ， 利 用 区 域 属 性 编码 来 表示 面 区 域 ， 
庄 廷 建立 了 距离 一 余弦 二 维 网 格 ， 通 过 获取 网 格 中 的 随机 点 数 ”通过 比较 区 域 属性 编码 来 度量 模型 之 间 的 相似 性 四。Jin 通过 推 
量 , 来 构建 三 维 模型 的 距离 -余弦 联合 形状 分 布 矩阵 。 利 用 分 布 。 导 拉 伸 形 变 能 量 的 三 维 版 本 ， 结 合 固定 边界 条 件 ， 将 表面 参数 
和 矩阵 之 间 的 LL; 距离 来 表示 模型 之 间 的 相似 度 , 并 采用 二 进 制 粒 比 扩展 到 体积 参数 化 ， 使 其 具有 拉 伸 最 小 化 的 优点 uM。 皇甫 中 
子 群 优化 算法 来 计算 模型 之 间 的 相似 性 外 。 张 开 兴 将 两 个 自 Sd nal nme 
面 的 相似 性 评价 问题 转换 为 两 个 曲面 局 部 特征 之 间 的 相似 性 ”来 表示 三 维 CAD 模型 , 使 用 图 谱 来 描述 局 部 特征 。 同 时 ， 

评价 问题 。 利 用 自由 曲面 参数 域内 生成 的 曲率 云图 来 提取 类 a 
Sift 算 子 ， 通 过 比较 局 部 特征 的 相似 程度 来 评价 自由 曲面 之 间 息 ， 生 成 多 个 零 忆 格 四 边 子 域 ， 并 使 用 边界 插值 算法 和 内 点 优 
的 相似 性 外 。 白 静 将 各 类 非 线 性 特征 统一 表示 为 属性 区 化 算法 来 实现 子 域 参数 化 03。 潘 万 彬 以 区 域 之 间 的 对 应 面 和 三 
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录用 稿 高 雪 瑶 ， 等 : Da (A 

维 尺寸 约束 图 为 基础 ， 将 简单 局 部 区 域 的 形状 和 定位 信息 传播 - Wi 

到 可 重用 区 域 131。 ul S(u,vi) ... 了 .9COyn) 
本 文 利 用 面 组 成 边 数 差 异 ， 来 构建 源 模型 与 目标 模型 之 间 RE 

的 面相 似 度 和 矩阵。 以 面相 似 度 矩 阵 为 基础 ， 使 用 粒子 群 算法 来 丰 SW) 

搜索 源 模型 与 目标 模型 之 间 的 最 优 面 匹配 序列 。 以 此 为 基础 ， a 

计算 源 模型 和 目标 模型 之 间 的 相似 性 。 um Snv) oo Sm bn) 


其 中 : 源 模型 的 面 序列 号 w1、u2、...、um， 作 为 面相 似 度 矩 阵 


1 ” 源 模型 与 目标 模型 之 间 的 面相 似 性 计算 


的 行 标 ， 目 标 模型 的 面 序列 号 vv、v2、...、vn， 作 为 面相 似 度 矩 
模型 是 由 多 个 面 组 成 的 。 组 成 面 的 形状 差异 ， 造 成 了 模型 。 阵 的 列 标 。 在 面相 似 度 矩 阵 SAs 中 ， 源 模型 面 wi 和 目标 模型 本 

ly 状 之 间 千 差 万 别 , 在 计算 源 模型 与 目标 模型 之 间 的 相似 性 时 ， ”vw 之 间 的 相似 度 值 为 SQi, vj)。 

需要 累积 这 两 个 模型 面 之 间 的 相似 性 。 如 果 模 型 之 间 面 的 组 成 源 模 型 面 与 目标 模型 面 之 间 的 相似 程度 影响 着 两 个 模型 

边 数 差距 很 小 ， 那 么 这 两 个 模型 的 相似 度 就 会 很 高 ， 如 果 模 型 。 间 的 相似 性 。 对 于 两 个 模型 而 言 ， 如 果 它 们 的 组 成 面 越 相 似 ， 


之 间 面 的 组 成 边 数 差异 很 大 ， 那 么 这 两 个 模型 之 间 的 相似 度 就 ” 那么 这 两 个 模型 之 间 的 相似 度 也 就 越 高 。 如 果 它 们 的 组 成 面 差 


很 低 。 源 模型 A 和 目标 模型 B 如 图 1 所 示 。 异 越 大 ， 则 这 两 个 模型 之 间 的 相似 度 也 就 越 低 。 
根据 式 (1)， 构 造 源 模型 A 与 目标 模型 B 之 间 的 面相 似 度 
和 矩阵。 在 面相 似 度 矩阵 中 ， 用 源 模型 A 的 面 序号 WW、w2、.…、 
7 表示 行 号 , 用 目标 模型 B 的 面 序号 yi、 v2、...、 WW 表示 列 号 。 
源 模型 A 和 目标 模型 B 的 面相 似 度 矩 阵 如 下 所 示 : 


V1 V2 V3 V4 V5 v6 v7 

wm 1 1 1 1 1 0.8 0.8 

WW W U3 uw 1 1 1 1 1 08 0.8 
wa 1 1 1 1 1 0.8 0.8 

源 模 型 A 目标 模型 B um 1 1 1 1 1 08 08 
Ss us 1 1 1 1 1 0.8 08 

图 1 源 模型 A 和 目标 模型 B us 0.8 08 08 08 08 1 1 

lu 0.8 08 08 08 08 1 1| 


源 模型 A 包括 七 个 面 : 面 w、w2、w3、u4、uUs、us 和 ur7。 面 
U1 与 面 w、us、us 和 wz 相 邻 接 ; 面 w 与 面 w、w3、ue 和 wy 相 在 源 模型 A 中 ， 面 ui 与 面 w2、us、ue 和 wz 相 邻接 ， 在 目 
邻接 ; 面 ww 与 面 w2、u4、us 和 wz 相 令 接 ; 面 u4 与 面 w3、Us、u6 标 模型 B 中 ， 面 vi 与 面 vw、vs、v6e 和 vw 相 邻 接 。 由 式 (1) 可 知 ， 
和 wz 相 邻 接 ; 面 us 与 面 wi、 ua、 ue 和 wi 相 邻 接 ; 面 wk 与 面 四、 在 源 模型 A 和 目标 模型 B 的 面相 似 度 和 矩阵 中 ， 第 wi 行 第 vi 列 
wu2、U3、U4 和 us 相 邻 接 ; 面 ww 与 面 wi、w2、u3、u4 和 ws 相 邻 接 。 的 值 为 1。 在 目标 模型 B 中 ， 面 ww 与 面 V、v2、v3、v4 和 vs 相 
面 ww 有 四 个 相 邻 面 ， 因 此 ， 面 ui 的 边 数 为 4。 同 理 ， 可 以 计算 邻接。 由 式 (1) 可 知 ， 在 源 模型 A 和 目标 模型 B 的 面相 似 度 算 
出 其 他 面 的 边 数 。 面 wz、w、w4 和 us 分 别 有 四 条 边 , 面 we 和 wz ” 阵 中 ,第 ui 行 第 w 列 的 值 为 0.8。 

分 别 有 五 条 边 。 


2 ”基于 粒子 群 算法 的 模型 面 匹配 


目标 模型 B 包括 七 个 面 : 面 V、 ww、V3、vV4、V、v6 和 v7。 
其 中 ， 面 vy、w、V3、w4 和 vs 是 四 边 形 ， 分 别 有 四 条 边 ; 面 ve 如 果 两 个 模型 的 组 成 面 越 相似 ， 则 这 两 个 模型 之 间 的 相似 
和 vw 是 五 边 形 ， 分 别 有 五 条 边 。 度 也 就 越 高 。 对 于 两 个 模型 而 言 ， 其 组 成 面 之 间 存 在 着 多 种 对 
利用 式 (1) 来 计算 源 模型 面 wu; 和 目标 模型 面 w 之 间 的 相似 应 方案 。 本 文 以 面相 似 度 和 矩阵 为 基础 ， 使 用 粒子 群 算法 来 搜索 
度 SC vw)。 源 模型 面 与 目标 模型 面 之 间 的 最 优 匹配 序列 。 根 据 最 优 面 匹 配 

| [num(u) — numv,)| 上 序列 ， 来 计算 两 个 模型 之 间 的 相似 性 。 

max(nundu),numv,)) 在 使 用 粒子 群 算法 搜索 面相 似 度 矩阵 的 过 程 中 ， 设 定 一 个 
其 中 : num 几 表示 模型 面 1 所 包含 的 边 数 ，max(x, 表示 x 与 yy 4d 维 的 搜索 空间 。 种 群 的 大 小 为 p; xi 表示 第 t 个 粒子 的 位 置 ， 
之 间 的 最 大 值 。 由 公式 (1) 可 以 看 出 : 若 两 个 模型 面 所 含 的 边 数 即 在 面相 似 度 矩阵 中 搜索 到 的 列 序号 组 成 的 向 量 ;vi 表示 第 1 个 
差异 越 小 ， 则 这 两 个 面 的 形状 就 越 相 似 。 使 用 SGwi, vj) 来 构造 源 粒子 的 飞行 速度 ， 在 第 于 次 迭代 过 程 中 ， 使 用 x 表示 第 


模型 A 与 目标 模型 B 之 间 的 面相 似 度 矩 了 泗 Sap， 其 中 : m 表示 
源 模型 面 数 ，n 表示 目标 模型 面 数 。 如 果 m 大 于 n， 则 将 源 模 ” 个 粒子 的 个 体 最 优 位 置 ， 对 应 的 适应 度 函 数值 为 feveswt; 利用 
型 与 目标 模型 互 换 。 即 始终 保持 面相 似 度 矩 阵 的 行 数 要 小 于 等 
于 列 数 。 


Xe si 表示 全 局 最 优 位 置 ， 对 应 的 适应 度 函数 值 为 axev。 利 用 
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Xpperr 4 和 Xk, 来 更 新 X11 和 速度 V4*" ， 其 计算 过 程 如 式 (2) 
和 (3) 所 示 。 


k+l 
Vi 


kk k Lk k ok 
= WV, +CN(Xppess — Xi ) + Cn (XGpess — Xr ) 2) 


X =X 十 V G3) 


其 中 :表示 当前 的 迭代 次 数 ; w 表示 惯性 权重 ; 
因子 ; rr、% 表 示 取 值 在 [0, 1] 内 的 随机 数 。 

在 每 次 迭代 过 程 中 ， 为 了 使 粒子 向 最 好 的 位 置 进行 偏 移 ， 
为 每 个 粒子 设 定 了 对 应 的 适应 度 值 。 第 t 个 粒子 在 面相 似 度 秆 
阵 SaB 中 搜索 到 的 位 置 向 量 x00),7(2),…,j(m))。 其 中 ，j() 表 
示 与 源 模型 A 的 第 i 个 面 所 匹配 的 目标 面 号 (二 1,2,.…, m)。 第 t 
个 粒子 的 适应 度 函数 的 计算 过 程 如 式 (4) 所 示 。 


cl、c2 为 学 习 


f= YsLi, 0)] 四 


基于 粒子 群 算 法 的 面 匹配 过 程 如 下 所 示 ; 

a) 根 据 式 (1) 构 造 源 模型 A 和 目标 模型 B 之 间 的 
矩阵 SAB。 

b) 初 始 化 种 群 大 小 p， 惯 性 权重 w， 和 迭代 次 数 s， 学 习 因 子 
初始 化 粒子 的 位 置 x 和 速度 w ( 关 1, 2, … n); 初始 化 当 
前 迭代 次 数 厂 0。 

o) 根 据 式 (4) 计 算 第 t 个 粒子 的 适应 度 函 数值 fAD)(=1, 2,.…， 


面相 似 度 


Cl 2C23 


n)。 


四 判断 2s， 如 果 为 真 则 算法 结束 ,输出 x5, ; 否则 执行 


eo) 如 果 1D<f ， 则 =D)，Xbpow =X 。 


f) 如 果 上 .<f, 则 人 =/ ， Ws RP past » 


g) 根 据 式 (2)(3) 来 更 新 第 t 个 粒子 的 速度 VV 和 位 置 xX" 。 


hb) 太太 1， 执 行 步骤 四。 
在 模型 面相 似 度 和 矩阵 Sap 中 ， 利 用 粒子 群 算法 来 搜索 两 个 
模型 之 间 的 最 优 面 匹配 序列 。 最 优 位 置 解 向 量 为 0(1), j(2),.…， 
Jj(m))。 可 以 知道 : 源 模型 A 与 目标 模型 B 的 最 优 面 匹配 序列 为 
(1,7(01))、 G2,7(2)、 .~ Cm, 7))。 
根据 最 优 面 匹配 序列 , 从 源 模型 A 与 目标 模型 B 的 面相 似 
度 和 矩阵 Sas 中 提取 对 应 的 元 素 。 通 过 累积 源 模型 面 和 目标 模型 
和 之 间 的 相似 性 ， 就 可 以 得 到 A、B 两 个 模型 的 整体 相似 性 
SmodlA.B)。 模 型 相似 性 Siwoaca,B) 的 计算 过 程 如 式 (5) 所 示 。 


wa Sl[i, 二 . ] 
这 i, j(i) ee 


min(m,n) 


匹配 序列 中 的 第 之 个 值 ; SI 为 画 


SS mod(A, B) 一 


相 


其 中 :(i, jQ)) 表 示 最 优 理 


似 度 和 矩阵 Sas 的 第 i 行 第 j() 列 的 数值 ，min(m, 站 为 m,n 的 最 小 
值 。 


3 ”实验 


本 文 在 处 理 器 为 5 的 Windows 7 操作 系统 下 , 使 用 Matlab 
来 编写 粒子 群 算法 ， 实 现 对 面相 似 度 矩 阵 的 搜索 。 在 实验 中 ， 
共 选 取 三 个 CAD 模型 。 目 标 模型 和 源 模型 如 图 2 所 示 。 目 标 
模型 是 四 棱锥 ; 源 模型 A 是 四 棱锥 ; 源 模型 B 是 三 棱柱 ; 源 模 


型 C 是 带 有 通 孔 的 圆 台 。 


.SO 
和 ~、 a 


Va 


目标 模型 源 模型 A 源 模型 B 。。 源 模 型 C 
图 2 目标 模型 和 源 模 型 
在 度量 源 模 型 与 目标 模型 之 间 的 相似 性 时 ， 一 般 都 采用 贪 


心 策略 来 寻找 这 两 个 模型 之 间 的 最 优 面 匹配 序列 。 为 了 验证 本 
文 所 提出 方法 的 有 效 性 ， 共 进行 了 两 组 对 比 实验 。 

在 实验 中 ， 利 用 式 (1) 来 构造 源 模型 与 目标 模型 的 面相 似 度 
矩阵 SAB。 分 别 使 用 贪心 算法 和 粒子 群 算法 来 搜索 这 两 个 模型 


之 间 的 最 优 面 匹配 序列 。 根 据 最 优 面 匹配 序列 ， 从 和 矩阵 SAs 中 
提取 源 模型 面 与 目标 模型 面 之 间 的 相似 性 ， 并 使 用 式 (5) 来 计算 


模型 之 间 的 相似 性 。 
标 模型 包括 五 个 面 : 面 vi、 面 vi 与 面 
VW、v4 和 vs 相 令 接 ; 面 :与 面 wm、 和 vs 相 邻 接 ; 面 vi 与 面 
多 、vV4 和 vs 相 邻 接 ; 面 ww 与 面 wm 、 和 vs 相 令 接 ; 面 vs 与 面 
7 yw 、 和 由 相 邻 接 。 面 有 三 个 相 邻 面 ， 因 此 ， 面 vi 的 边 
数 为 3。 同 理 ， 可 以 计算 出 其 他 面 的 边 数 。 面 w、v3 和 va 分 别 
有 三 条 边 ， 面 vs 有 四 条 边 。 
以 图 2 中 的 源 模型 A 为 例 , 使 用 粒子 群 算法 来 搜寻 源 模型 
A 和 目标 模型 之 间 的 最 优 面 匹配 序列 。 源 模型 A 包括 五 个 面 : 
其 中 ， 面 ww、u2、w3 和 ws 为 三 角 
分 别 有 3 条 边 ; 面 us 是 四 边 形 ， 含 有 四 条 边 。 
使 用 式 (1) 计 算 源 模型 面 与 目标 模型 面 之 间 的 相似 性 ， 并 构 
造 源 模型 A 与 目标 模型 之 间 的 面相 似 度 和 矩阵 ， 如 下 所 示 : 
V1 v2 V3 V4 V5 ] 
Ul 1 1 1 1 0.75 
uz 1 1 1 1 0.75 
U3 1 1 1 1 0.75 


u4 1 1 1 1 0.75 
us 0.75 0.75 0.75 075 1 | 


在 源 模 型 A 中 ， 面 ui 与 面 w2、u4 和 us 相 邻 接 ; 在 目标 模 
型 中 ， 面 vi 与 面 y2、v4 和 vs 相 邻 接 。 由 式 (1) 可 知 ,在 源 模型 A 
和 目标 模型 的 面相 似 度 和 矩阵 中 ， 第 i 行 第 vi 列 的 值 为 1。 在 目 
标 模 型 中 , 面 vs 与 面 Vi、v2、v3 和 va 相 邻 接 。 由 式 (1) 可 知 ,在 源 
模型 A 和 目标 模型 的 面相 似 度 矩 阵 中 ， 第 ui 行 第 vs 列 的 值 为 


ww、 VW、V4 和 vs。 
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面 wi、u2、 U3、uU4 和 us。 
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以 源 模型 A 和 目标 模型 的 面相 似 度 矩阵 为 基础 ， 使 用 粒子 
群 算法 来 搜索 源 模型 A 与 目标 模型 之 间 的 最 优 面 匹配 序列 。 得 
到 最 优 面 匹配 序列 为 : (1, 3)、(2, 2)、(3, 1)、(4, 4)、(5, 5)。 同 
时 ， 使 用 式 (5) 来 计算 源 模型 A 与 目标 模型 之 间 的 整体 相似 性 ， 
其 数值 为 1。 
源 模型 A、 源 模型 B 和 源 模 型 C 与 目标 模型 之 间 的 相似 性 
计算 结果 如 表 1 所 示 。 
表 1 源 模型 与 目标 模型 之 间 的 相似 性 
源 模型 ” 面 数 ” 边 数 ” 顶 点数 ”贪心 算法 ”粒子 群 算法 


让 
十 


A 5 8 5 1 Eb 
B 5 9 6 0.9 0.9 
C 4 8 0 0.75 0.8125 


从 形状 上 看 , 源 模型 A 与 目标 模型 完全 一 致 ， 其 相似 性 数 


值 为 1， 源 模型 B 与 目标 模型 的 相似 度 较 大 ， 其 相似 性 数值 为 
0.9; 源 模型 C 与 目标 模型 有 一 定 的 相似 性 ， 其 相似 性 数值 为 
0.812 5 。 


从 表 1 可 知 , 源 模型 A 和 目标 模型 是 同一 个 模型 ， 贪 心算 
法 和 粒子 群 算法 获得 的 数据 是 一 致 的 。 对 于 源 模型 B， 贪 心算 
法 和 粒子 群 算法 获得 的 数据 也 是 一 致 的 。 对 于 源 模 型 C， 它 是 
一 个 关键 性 模型 。 使 用 粒子 群 算法 所 获得 的 相似 性 数值 比 贪心 
算法 要 高 ， 其 结果 提高 了 8.33%。 从 形状 上 看 ， 源 模型 C 与 目 
标 模 型 有 一 定 的 相似 性 。 相 对 于 贪心 算法 而 言 ， 本 文 所 提出 的 
方法 更 能 准确 地 度量 两 个 模型 之 间 的 相似 性 。 
在 实验 中 ， 使 用 贪心 算法 搜索 面相 似 度 和 矩阵 ， 其 时 间 复 杂 
度 为 O(mz)。 使 用 粒子 群 算法 搜索 面相 似 度 矩阵 ， 其 时 间 复 杂 
度 为 Opsmn)。 尽 管 粒子 群 的 时 间 复 杂 度 要 高 于 贪心 算法 ， 但 
是 粒子 群 算法 在 度量 模型 之 间 的 相似 性 时 ， 效 果 要 好 于 贪心 算 


4 ”结束 语 


本 文 利用 面 的 边 数 差异 来 构建 面相 似 度 矩阵 。 通 过 
算法 来 搜索 源 模型 与 目标 模型 之 间 的 最 优 面 匹配 序列 。 以 最 优 
面 匹配 序列 为 基础 ， 通 过 累积 源 模型 面 与 目标 模型 面 之 间 的 相 
似 性 来 度量 两 个 模型 的 整体 相似 度 。 实 验 粒子 群 算 


粒子 群 


吉 果 表明 : 


度量 模型 之 间 的 形状 差异 。 
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